function F      =   FIGs_C1_C2_C3_C4_find_peg_shocks(guess)

global y0 long drop rhoA  replic iorder prune p_int sug peg dr

e1              =   [0 0 0.01 guess(1) 0 0 guess(2) guess(3) guess(4)...
                     guess(5) guess(6) guess(7) guess(8)]';

irf_c           =   FIGs_C1_C2_C3_C4_irf_dynare(dr, e1, long+1, drop, replic, iorder);

F               =   irf_c(p_int,1:peg)